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Abstract 

We develop a relativistic Langevin dynamics under the background of strongly interacting quark- 
gluon fluid described by the (3+l)-dimensional hydrodynamics. The drag force acting on charm 
and bottom quarks is parametrized according to the formula obtained from the anti-de-Sitter 
space/conformal field theory (AdS/CFT) correspondence. In this setup, we calculate the nuclear 
modification factor Raa for the single-electrons from the charm and bottom quarks to extract the 
magnitude of the drag force from the PHENIX and STAR data. The Raa for single-electrons 
with pr > 3 GeV indicates that the drag force is close to the AdS/CFT prediction. Effects of 
the drag force on the elliptic flow of single-electrons are also discussed. Moreover, we predict 
the electron-muon correlation which is closely related to the heavy-quark pair correlation in hot 
matter. 



1. Introduction 

The properties of quark-gluon plasma (QGP) are actively investigated by relativistic heavy 
ion collision experiments at Relativistic Heavy Ion Collider (RHIC) and will be studied at the 
Large Hadron Collider (LHC) in the future. The space-time evolution of the created hot medium 
is described nicely by the relativistic hydrodynamic model in which vanishing viscosity is as- 
sumed. This fact implies the strongly interacting nature of the QGP even above the deconfine- 
ment transition temperature T c ~ 170 MeV. In heavy ion collisions, majority of the particles, 
namely light soft particles such as up, down and strange quarks and gluons, constitute the hot 
medium whose space-time evolution is described by the hydrodynamic model. However, heavy 
quarks, namely charm and bottom quarks, do not join the medium but act as impurities in the 
medium because they have much heavier masses (m c = 1.5 GeV, m/, = 4.8 GeV) than typical 
temperature available at RHIC (T ~ 300 MeV). Not only heavy quarks but also jets and J/^s 
act as impurities in the medium. 

In this report, we focus on heavy quarks in the QGP, especially on their diffusive dynamics 
in the medium and on its influence on the spectra of their observables in experiments, namely 
single-leptons decayed from heavy mesons. We perform simulations of heavy quark diffusion |lj] 
in the background of expanding hot medium and of the subsequent heavy quark hadronization 
followed by semileptonic decay and make a quantitative comparison with experimental data H. 
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2. Langevin + Hydrodynamics Model 



Since masses of heavy quarks are much larger than the typical temperature of surrounding 
medium while kinematics of heavy quarks can be relativistic in actual experiments, the heavy 
quark dynamics may be described by the relativistic Langevin equation. The basic equations in 
the rest frame of matter are 

Ax = -At, Ap = -T(p)pAt + |(0, (1) 
E 

<6(0£/(O> = D ij {p)5 tt ,At, (2) 

where we use the specific parametrization, Tip) - yT 2 /M and Dy(p) - 2yT 3 (E + T)6ij/M, 
motivated by the drag force calculated on the basis of the anti-de-Sitter space/conformal field 
theory (AdS/CFT) correspondence [3]. According to an attempt to translate the drag parameter 
y of the supersymmetric Yang-Mills plasma into the one in strongly interacting QGP, the drag 
parameter is estimated to be y — 2.1 + 0.5 [|4J]. 

The combination of the relativistic Langevin dynamics of heavy quarks and hydrodynamic 
expansion of hot medium is implemented by the following procedures, (i) We first specify the 
initial condition of heavy quarks. The initial distribution is given by multiplying the nucleon- 
nucleon collision by the number of binary collisions, namely the Glauber model, (ii) Heavy 
quarks diffuse in the expanding medium. The local flow vector and temperature are given by the 
hydrodynamic model, (iii) Heavy quarks hadronize to heavy mesons, namely D or B mesons, via 
independent fragmentation when confinement phase transition takes place, (iv) Heavy mesons 
decay to single-electrons via semileptonic decay. 

Some comments are in order here: (a) The initial heavy quark momentum distribution is 
computed by PYTHIA event generator, which is based on leading order perturbative QCD and 
is found to be reliable only at high transverse momentum pj region. For this reason, we limit 
ourselves to the comparison of the single-electron results from our simulation with those from 
experiment only in their spectral shape at p? > 3.0 GeV. (b) Since we assume first order con- 
finement phase transition in the hydrodynamic model, the mixed phase of QGP and hadronic 
phase causes ambiguity in the determination of the time when heavy quarks hadronize to heavy 
mesons. We treat this ambiguity as a systematic error by performing simulations with different 
freezeout criteria. 



3. Raa for single electrons and electron-muon azimutal correlation 

We show in FigQ]the nuclear modification factor Raa of single-electrons defined as, 

1 dN'JdpT 

R aa (pt)=— * +A / ■ (3) 
A'con dN* +v /dpT 

The collision geometry for Raa is fixed by impact parameter b = 3.1 fm and the rapidity range 
considered is \t] p \ < 0.35. Numerical results of simulation with different drag parameter y and 
with different choices of freezeout condition are plotted. By comparing our simulation with 
experimental data at pj > 3.0 GeV, we find that the drag parameter y — 1-3 is favored. Note that 
this value is consistent with the prediction from the AdS/CFT correspondence y = 2.1 + 0.5. 

The degree of thermalization of heavy quarks can be estimated by the ratio of two time scales, 
the stay time t$ and the relaxation time tq. With a drag parameter y = 1-3, ts is 3-4 fm and tq 
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is 2-7 fm for charm and 7-20 fm for bottom quarks [ 1 ] . Therefore the charm quarks are partially 
thermalized while bottom quarks are not thermalized at RHIC. 

As for elliptic flow V2 of single-electrons, our simulation and experimental data still suffer 
from poor statistics at pj > 3.0 GeV and definite conclusions cannot be drawn at the moment. 
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Figure 1 : Comparison of the nuclear modification factor i?AA of single-electrons in our hydro + heavy-quark model and 
the experimental data. The collision geometry is fixed by the impact parameter 3.1 fm and the single-electrons with mid- 
rapidity \r/ p \ < 0.35 are counted. Results of our simulation with different drag parameters and with different freezeout 
conditions are drawn. 

We show in Fig|2]our results of electron-muon azimuthal correlation. This azimuthal correla- 
tion is a clean observable and free of the contamination by di-leptons from vector meson decays. 
The trigger particles are single-electrons with pj > 3.0 GeV and the associated particles are 
single-muons with pj > 3.0 GeV. The high pj trigger electron is desirable because it takes over 
the direction of its parent heavy meson with even higher pr, which also inherits the direction 
of the heavy quark just before hadronization. Thus the azimuthal correlation between a heavy 
quark and an anti-heavy quark pair at their freezeout can be studied by the correlation between 
a single-electron and a single-muon. Collision geometry is again fixed by the impact parameter 
b = 3.1 fm. We show our results with the pseudorapidity region \rj\ < 1 .0 in Fig|2] In this corre- 
lation, we can see only one peak in the away side because a single-electron and a single-muon 
cannot be produced simultaneously in one semileptonic decay. For the drag parameter y > 3, the 
away side quenching is significant. Therefore, clean measurement of the azimuthal correlation 
between a single-electron and a single-muon with high transverse momentum pj > 3.0 GeV can 
probe the strongly interacting QGP quite well. 



4. Summary 

In this report, we showed our studies on the heavy quark diffusion in the dynamical QGP fluid 
on the basis of the relativistic Langevin equation combined with the relativistic hydrodynamics. 
Motivated by the formula from the AdS/CFT correspondence for strongly interacting plasma, 
we parametrized the drag coefficient by T = yT 2 /M with the dimensionless coefficient y as a 
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Figure 2: Azimuthal correlation of a single-electron and a single-muon decayed from a heavy quark and anti-heavy 
quark pair. Collision geometry is fixed by the impact parameter 3.1 fm. Both trigger electrons and associate muons have 
transverse momentum pj > 3.0 GeV and are in mid-pseudorapidity \rj\ < 1.0. 

parameter to be extracted from experiment. We obtained single-electron spectra by solving the 
Langevin equation for heavy quarks on the space-time evolution of background fluid and by 
simulating their hadronization and subsequent semileptonic decay. 

The comparison of /?aa for pj > 3 GeV suggests that the drag coefficient could be as large 
as y — 1.0 - 3.0, consistent with the prediction by the AdS/CFT. We have also found that az- 
imuthal electron-muon correlation with high transverse momentum pj > 3 GeV has a significant 
dependence on the drag parameter and thus can serve as a good probe for the strongly interacting 
QGP. 
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